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We present a phenomenological calculation of the quasiparticle-interference (QPI) pattern and inelastic Neu- 
tron scattering (INS) spectra in iron-pnictide and layered iron-selenide compounds by using materials specific 
band- structure and superconducting (SC) gap properties. As both the QPI and the INS spectra arise due to 
scattering of the Bogolyubov quasiaprticles, they exhibit an one-to-one correspondence of the scattering vectors 
and the energy scales. We show that these two spectroscopies complement each other in such a way that a com- 
parative study allows one to extract the quantitative and unambiguous information about the underlying pairing 
structure and the phase of the SC gap. Due to the nodeless and isotropic nature of the SC gaps, both the QPI 
and INS maps are concentrated at only two energies in pnictide (two SC gaps) and one energy in iron-selenide, 
while the associated scattering vectors q for scattering of sign-changing and same-sign of the SC gaps change 
between these spectroscopies. The results presented, particularly for newly iron-selenide compounds, can be 
used to test the nodeless d— wave pairing in this class of high temperature superconductors. 
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I. INTRODUCTION 

The most crucial information to unravel the mechanism of 
pairing is the structure of the superconducting (SC) gap func- 
tion, a measure of the amplitude and the phase of electron 
pairs. Although the SC-gap function of conventional phonon- 
mediated (attractive pairing interaction) superconductors has 
the same sign all over the /c— space (s— wave symmetry), 
that of spin-fluctuation mediated (repulsive pairing interac- 
tion) superconductors, within the weak to intermediate cou- 
pling scenario,^ is expected to exhibit a sign reversal between 
the Fermi momenta connected by the characteristic 'hot-spot' 
wave vector Q of spin-fluctuations. ^^1 As a consequence of 
the sign reversal, nodal planes in which the SC gap van- 
ishes should exist in k space. If such nodal planes inter- 
sect the Fermi surface (FS), as in (i— wave cuprate, nodal 
quasiparticle states emerge in the low-energy spectra. How- 
ever, if the FS consists of small pockets and do not touch 
the nodal line, as pairing in iron-pnictide and single 
layer iron-chalogenides"^ or d-wave pairing in double layered 
iron-selenide (Fe2Se2),^ ^ develop despite a sign reversal of 
the pairing. Therefore, to precisely establish an unconven- 
tional pairing symmetry and separate it from a conventional 
5— wave, the relative sign of the SC gap on the FS must be 
established. 

The magnetic resonance mode develops in the SC state 
and was shown theoretically to arise from the sign-flip of the 
SC gap at the 'hot-spot' P^Jli^ However, there are other the- 
ories which also reproduce the resonance mode without the 
requirement of sign change of the SC gap.fJ^Jf^ Alternative 
and important new spectroscopy comes from the quasiparti- 
cle interference (QPI) pattern, measured by scanning tunnel- 
ing microscopy (STM), which in principle visualizes all pos- 
sible elastic scattering vectors. It therefore can distinguish 
the sign-changing 'hot-spot' vectors by studying its evolution 
with varying magnetic Magnetic field breaks the 

time-reversal symmetry of the SC quasiparticles and thus il- 



luminate those scattering vectors which scatter quasiparticles 
of same pairing phase, thanks to the remarkable properties 
of the Bogolyubov coherence f actors. These observations 
can be questioned based on at least two arguments: (1) mag- 
netic field induced vortex states drastically redistribute the 
spectral weight across the 'bright- spots' on the constant en- 
ergy surface, making this procedure complicated; (2) a mag- 
netic impurity always carries a scalar potential which allows 
simul taneo us scattering of quasiparticles of opposite pairing 
phase.ES^We show in this paper that the INS and QPI maps 
complement each other and allow determination of the rel- 
ative phase. Thus, we propose, a comparative study between 
the Inelastic neutron Scattering (INS) spectra and the QPI pat- 
tern, taken together, can be a viable tool to quantitatively and 
unambiguously determine the relative sign of the SC gap. 

Our approach is based on a very simple observation: Both 
the INS and the QPI patterns arise from the similar scattering 
of the Bogoliubov quasiparticles, forming Cooper pairs (in- 
elastic and elastic scatterings, respectively) and thus their ob- 
served dispersions, l^(q)s, must have an one-to-one relation- 
ship in these two spectroscopies. In INS spectra, a scattering q 
vector will be observed if it connects two quasiparticle states 
at and kj at which sgn [ Ak^ ] ^ sgn [ Ak ^ ] and thus its 
energy scale will be determined by | Ak^ | + | Ak^ | . The same 
vector will only appear in the QPI spectra if | Ak^ | = | Ak^ |, 
because QPI probes elastic scattering. The comparison be- 
tween the INS and the QPI maps can therefore emphasize the 
regions of momentum space with opposite sign of the gap 
function, connected by scattering momenta. This emphasis 
can be further facilitated by studying the magnetic field de- 
pendence of the QPI maps. Such a comparative study has 
proven to give valuable information about the pairing symme- 
try in cuprates.'^ 

Iron-pnictide:- In the case of iron-pnictide superconduc- 
tors, however, the presence of multiple bands at the Fermi 
level (Ef) and multiple SC gap amplitudes at each bands 
make the aforementioned analysis much more complicated 
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and exotic. In these compounds, the FS consists of discon- 
nected two concentric hole pockets (namely, a and {3 pock- 
ets) and one electron pocket (7 pocket) centered at T and 
M points, respectively. Prominent nesting between the hole 
pockets and the electron pocket has been shown theoretically 
to lead to a s^— pairing in this family of superconductors .El 
Within this pairing, all the FS pockets possess nodeless and 
isotropic SC gap, but the sign of the gaps is reversed between 
the hole pockets and the electron pocket. To complicate the 
story, the tw o hole p ockets acquire very different magnitude of 
the SC gapsi ^^ l ^^ l S l Therefore, macroscopic geometrical sep- 
aration of the two phases, as in the cuprate corner-junction 
experiment, can not be performed in these systems to test the 
pairing symmetry.l^ 

In our earlier study, we have theoretically shown that be- 
cause of the two values of the SC gap in two hole-pockets, 
the INS spectra is split into two energy scales. Furthermore, 
as the two hole pockets have very different area as a function 
of doping, the corresponding q vectors are also different at 
these two resonances. Here, we show that these INS scatter- 
ing vectors also appear in two different energy scales of the 
QPI patterns. Additional scattering vectors between the states 
of same sign of SC gaps, which were prohibited in INS, ap- 
pear in QPI. These results are consistent with STM data on 
electron doped Ba(Fei_a;COa;)2As2 as a function of doping.Sl 
Furthermore, the evolution of the QPI pattern as a function of 
magnetic field is also consistent with STM studies in single 
layer Fe(Se,Te) superconductors.^^ 

Fe2Se2 superconductors:- The recent discovery of high- 
Tc superconductivity in double layered Fe2Se2 based super- 
conductors makes the above story even more interesting and 
exotic. ^^ ^^ The FS in this family of superconductors only 
hosts electron pockets at M point (no hole-pocket is present 
here as in iron-pnictide discussed above). Therefore, the lead- 
ing nesting vector occurs between the two electron pockets 
which lead to a (i-wave pairing {dj.2_y2 in 1 Fe unit cell and 
dxy in 2 Fe unit cell notation).!^ Unlike the d-wave pairing in 
cuprates, here it gives rise to nodeless and isotropic SC gap 
on the FS, in consistent with experiment J^^lEl gap func- 
tion changes sign between the two electron pockets. We and 
others^ have shown earlier that such a d-wave pairing leads 
to a magnetic resonance mode at Q which is observed later by 
INS measurements.^^ 

Here we show that this INS scattering vector Q will also 
show up in the QPI vector at the bias energy equals to the 
SC gap magnitude. Furthermore, we phenomenologically 
demonstrate that with the application of the magnetic field, 
other QPI scattering vectors at which the SC gap do not 
change sign can be illuminated. The relative evolution of the 
QPI maps as a function of magnetic field will provide a valu- 
able test for the relative phase of the SC state on the FS in this 
family of superconductors. 

The paper is organized as follows. In Sec. II, we provide 
the formalism for both QPI and INS spectra. We extract out 
very simple equations for these two spectroscopies which can 
be used with experimental inputs to reconcile them. In Sec. 
Ill, we present the results for iron-pnictide, whereas the same 
results but for layered iron-selenide are given in Sec. IV. The 



magnetic field dependence of the QPI map is computed for 
layered iron-selenide system in Sec. IVA. Finally, we con- 
clude in Sec. V. 



II. FORMALISM 

A direct correlation between the INS and the QPI spectra 
can only be made within the bare level in which the same 
green's functions are involved. This correspondence there- 
fore would at least qualitatively be correct for dressed Greens 
functions. In fact, one can again evaluate the Green's func- 
tions for the filled state from ARPES spectral weight A, as 
G(k, iujn) = / A(k, iuo')/{iuJn — 00'), which is exact for a 
single band case, but averaged over orbital indices in a multi- 
band system. In a multiband superconductor, the Green's 
functions for the normal and anomalous part can be written 
in the eigenbasis as 

^ (icjn - Ej,{k) iujn + Ej,{k) ) ' 

respectively, n is the Matsubara frequency. Here Eiy{k) = 
±y^^^^J^)~^^~^Jk) is the Bogolyubov quasiparticle 
band, ^^(k) is the non-interaction band modeled by tight- 
binding parametrization to the material specific LDA bands, 
as discussed later. Aj,(k) = A^^(k) is the momentum de- 
pendent SC gap with is the band specific SC gap mag- 
nitude, ^(k) is the structure factor for the pairing symmetry 
which is taken to be same for all bands in a given system. 

a,(k)(/3,(k)) = are the Bogolyubov co- 

herence factors for the quasiparticle states ±Ejy{k), respec- 
tively. The above Green's functions can be projected to the 
orbital basis as Gsp(k,'^^n) = (0f,(k)|Gi.(k, icjn)|0i^(k)) 
and same for F. 0f,(k) represent the eigenvectors of band 
onto s^^ orbital. 

QPI profile: -STM measures local density of state (LDOS) 
of quasiparticles which arrive to the tip after going through 
multiple intrinsic elastic scattering (due to magnetic and non- 
magnetic scattereiE3I3^ in the system. Such scattering vec- 
tors can be visualize d by Fo urier transforming the LDOS into 
the q space to 

X \V^\[G^k,^n^)G^\k^cl,^nr^) 
- F^(k, inn)F^' (-k - q, -iQn)] . (3) 

Here M is the orbital to band matrix element made 
of the eigenvectors as M!^^^^{k,q) = ^I^*(k + 
q)^^/(k)^^*(k)^^(k + q). The matrix-element can be 
important for shaping the spectral weight distributions 
for multiband case. Nevertheless, the relative intensity 
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of the scattering vectors is determined by the nesting 
conditions and the SC gap ampHtude for elastic scatter- 
ing. In addition, a scattering matrix-element, C(k, q), 
appears due to the coherence factors, ^^^(k) and /3j,(k), 
of the Bogolyubov quasiparticles which is sensitive to the 
momentum-dependent phase of the SC order parameter 
and the symmetry of the scattering potential, Vq.^^^^ The 
green's function is expressed in terms of the spectral weight 
asG-(k,zl]n) = -J^^^A-{k,uj')/{inn-oj'). If no 
many-body broadening is included then the spectral weight 
can be presented by delta function. Taking these facts into 
account, we simplify Eq. [3] to calculate the total QPI map 
as,Ci 



5(q,(^QPi)^|Fc 



Cuu'iK q) 



xs{nQPi - E,{k))s{nQPi - E,{k + q)). 

(4) 

[Here we assumed M = 1 for simplicit>^and took analytical 
continuation i^ln ^qpi + '^V (V is a small broadening).] 
The explicit form of the scattering matrix element as given 
inRef.[l6]and[l7|isa^/(k,q) = (sgn[A^(k)]sgn[A^/(k + 

q)]a^(k)a^/(k + q) T /3^(k)/3^/(k + q))^ where the 
negative and positive signs represent scattering through 
a scalar (even under time-reversal) and a magnetic (odd 
under time-reversal) potential, respectively. Further- 
more, scattering from disorder that converts electrons into 
holes as they are scattered, gives rise to the coherence 
factor C^^/(k,q) = (sgn[A^(k)]a^(k)/3^/ (k + q) + 
sgn[A,. (k + q)]/3, (k)a,. (k + q)) (sgn[A,(k)]sgn[A,. (k + 
q)]a,(k)a,.(k + q) + /3,(k)/3,.(k + q)). At Ef the 
two delta function in Eq.|4] gives l]Qpj(q) = |Aj,(ki?)| = 
|Ajy/(ki? + q)|. Imposing this constraint, we obtain 
a^yikp) = a^/(kF + q), and /3^(kF) = f^u^ik^ + q) at 
Ep (by substituting = in the equations for a, /3 given 
above). This leads to a systematic selection rule for q vectors 
which depends on the nature of the scatterer: In the case of 
weak scalar potential scattering, C for those q at which 
sgn [Aj,(ki?)] = sgn [A^^/ (ki? + q)] . By contrast, for scatter- 
ing off magnetic impurities or gap inhomogeneities, C ^ 
for those q at which sgn [Ai^(ki?)] ^ sgn[Aj^/(ki7 + q)] . 
Therefore, focusing on the low-energy region where 
^QPi < A, we can write the conditions to obtain a non- 
vanishing QPI vector as 



which satisfies Eq. [6] as well but its intensity will depend on 
the relative strength of the potential. We take V^(q) as a con- 
stant potential. 

INS spectra: -The calculation of the INS spectra follows 
similarly with the exception that the latter is an inelastic scat- 
tering of the quasiparticle spectra. INS probes the imaginary 
part of the susceptibility which can be written in the orbital 
basis as^^ 



^^QPi(q) 



A,(kf)| = |A,,(kF + q)|, (5) 
sgn[Aiy(ki7)] 7^ sgn[A^/(ki? + q)] for scalar imp., (6) 
sgn[A^(kF)] = sgn[Ai,/(kF + q)] for mag. imp. (7) 



Note that as magnetic impurity is always associated with 
a scalar potential, it, in principle, involves QPI scatterings 



XOrstui^.^ 



X [ff'{k,iuJn)G''\k^q,iUn^inm) 
+ F^(k, iuJn)F^' (-k - q, -iun - 



(8) 



We have shown earlier in Ref. 12 that in many cases the dif- 
ference between realistic matrix element and matrix element 
assumed to be a smooth function of energy and momentum 
and thus is not important in calculating the magnetic scatter- 
ing structure. The INS dispersion is mostly governed by the 
locus of the discontinuous jumps in xo and due to Kramers- 
Kroning relationship, corresponding Xo attains a peak at the 
same location [within random-phase approximation (RPA), 
the peak position shifts slightly to a lower energy]. Similar 
to Eq. |4] above, we absorb the INS scattering matrix-element 
term in Cj^j^' (k, q) and performing the Matsubara summation 
in Eq.[8] we obtain the total ECS xo as 

X(5(aNs-^.(k)-^.(k + q)). (9) 

Here, the explicit form of C is Ciyjy/(k, q) = I3jy{k)aiy'{k + 
q)[a,(k)/3,.(k + q) - /3,(k)a,.(k + q)] = (l - 
sgn[A^(ki.)]sgn[A^.(k^ + q)]) at^;^. 

This implies that the magnetic structure in BCS Yn below 
^INS ^ 2 A is entirely governed by two conditions P'^'SII 



INS 



s(q) = |A^(kF)| + |A^/(kF+q)|, 



sgn[A^(ki.)] ^sgn[A^.(kF + q) 



(10) 
(11) 



We will use QPI Eqs.[5][7]and INS Eqs.[TOl&[TT]to perform 
a comparative analysis of the two spectroscopic data in iron- 
pnitide and iron-selenide superconductors. 



III. RESULTS ON IRON-PNICTIDE 

FSs in pnictide: -The low-energy Hamiltonian of Iron- 
pnictide system is dominated by five (i-orbitals of the Fe 
atoms. We take the tight-binding model from Ref. 38 where 
the parameters are fit to the corresponding first-principles dis- 
persion. The doping is evaluated within rigid-band shift ap- 
proximation. At X = —0.1, the computed FS consists of two 
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FIG. 1. (a) Computed FS for iron-pnictide within five band tight- 
binding model at a representative electron doping of x = 0.1. The 
BZ is chosen for 1 Fe unit cell notation. The black to white back- 
ground depicts the -pairing symmetry which takes the form of 
2 cos (kxa) cos (kya) in the 1 Fe unit cell. The arrows give different 
interband scattering channels which constitute QPI and INS maps, 
(bl) According to Eq. [s] the QPI map at |A/3| < |Aa| ^ |A^| 
only reveals the intraband scattering (schematic) within the /3— FS 
pocket. (b2) All possible interband scatterings from a 7-FS (qi) 
and from 7 7-FS, (q3,4), but not from /3 'j-FS (€[2) become 
turned on at Hqpi. (c1)-(c2) Due to inelastic scattering process in 
INS, the scattering between a- or /3- to 7-FSs is only allowed which 
changes sign of the SC gap, according to Eq. [TT] q2 for /3 7-FS 
appears at ^^Jns — 1^/3 1 + 1^7 1 in (cl). (c2) qi for a 7-FS 
appears at QiNS = |Ac| + |A^|. As |Ac| ~ |A^|,INS spectrum in 
(c2) will resemble QPI map in (b2) at zero magnetic field. Note that 
in the INS calculation, we have not included the Umklapp scattering 
which will symmetrize the INS spectra with respect to the Umklapp 
vector Q. This is done to facilitate the direct comparison with QPI 
maps which are not calculated using Umklapp scattering to mimic 
the experimental procedure. 



concentric hole pockets at F points which are called a- [in- 
ner pocket as depicted by red line in Fig. 1(a)] and /3-pockets 
(blue line) and one electron pocket at M point (green line). 
All results in the this paper are presented in the unfolded BZ 
coming from 1 Fe unit cell. 

SC gap properties .--From the shape of the FS pockets, there 
are at least four interpocket scattering channels exist in pnic- 
tide which span along various high-symmetry q directions 



as shown by arrows of different colors in Fig. 1(a). Among 
the four vectors, Mazin et al^ have shown theoretically that 
the nesting for qi and q2 are the strongest which lead to a 
sign-changing 5^ -pairing symmetry in this class of supercon- 
ductors. This phase symmetry of the SC gap is consistent 
with the spin- fluctuation mechanism of electron-pairing. The 
pairing and SC gap amplitudes have three essential properties 
which are relevant to our present study: (1) SC gap changes 
sign between the electron and hole pockets (black to white 
background colors in Fig. 1(a) reflect the 5^ -pairing symme- 
try), (2) SC gap magnitude on all FS pockets is nodeless and 
isotropic, and (3) evidence from ARPES^^, STM^^ and nu- 
merous bulk probes^^ indicate that. 



2\A, 



(12) 



at all dopings and for both electron and hole dopings. 

Sketch of QPI and INS maps:- All the aforementioned FS 
and SC gap properties lead to very different QPI and INS 
properti es in pnictide, than the ones obtained in single band 
cuprates.C^ 

(i) Nodeless and isotropic nature of the SC gaps make all 
the QPI and INS maps to concentrate at only two energy 
scales, instead of a characteristic dispersion seen in nodal and 
anisotropic d-wave gap in cuprates. QPI maps will be promi- 
nent only at ^qpi = |A/3| and ^qpi = |Ac,| ~ l^^], 
obeying Eq. p] While, the INS maps will show up only at 



Eq. 10 



1^/3 [+ ho 



and Qj^^ 



Aq, + , according to 



(ii) At ^^Qpi = |A/3|, no interband elastic scattering is 
allowed as |A/3| < |Aq,^^|. Therefore, the QPI map, in 
Fig. l(bl) only shows intra- /3-FS scattering which concen- 
trate near q = 0. As SC gap does not change sign on each 
pocket, some finite magnetic field will be necessary to illumi- 
nate these small q vectors. Of course, in real material, finite 
broadening can introduce some quasiparticle states of a, 7 
bands near I^qpi = | A/3 1 to visualize weak intensity at q2 (at 
zero magnetic field). 

(iii) At ^^Qpi = |Aq,| ^ |A^|, the interband scattering 
vector qi, intraband scattering vectors qa and q4 appear in 
Fig. I(b2). qs and q4 become illuminated at zero magnetic 
field while qi will gain more intensity at finite magnetic field, 
in consistent with the data from Ref. 17 

(iv) At ^^iNs = 1^/3 1 + I ^7 1' vector can be observed 
through INS study in Fig. l(cl). Note that, among the four 
maps shown in Fig. 1, this is the only place where q2 can be 
determined precisely. 

(v) Similarly, at ^^i^s = + I^tI' ^1 vector will be 
observed as sketched in Fig. I(c2). At zero magnetic field, 
this INS spectra will match exactly with QPI pattern shown in 
Fig. I(b2). 

(vi) Doping dependence of the q vectors and their energy 
scales (not studied here) can also be used to gain confidence 
on the FS topology and the location of the sign reversal of the 
pairing symmetry. The INS energy scales, ^^j^g obtain dome 
like behavior with doping in accord with the dome like behav- 
ior of the SC gaps as calculated in Ref. 12 By implication, 
the same doping dependence is expected in l^Qpj. Even in 
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a multi SC gaps pnictide system, such one-to-one correspon- 
dence is po ssible as all gaps show similar dome-like doping 
dependenceP^I^^ The area of the hole pockets increases with 
hole doping; simultaneously the same for the electron pocket 
decreases. This doping dependence is reversed for electron 
doping. Therefore, all q vectors will characteristically fol- 
low the doping dependence of their corresponding connecting 
FSs.fJ^ Interestingly, near x = 0.15 of electron doped side, 
a pocket disappears .^^^^ Therefore, qi should also disappear 
and hence ^^Qpi (at zero magnetic field), ^^i^g. These results 
are yet to be confirmed experimentally. 



A. Computed QPI maps of iron-pnictide 

We compute the QPI maps for pnictide using Eq.|4]for the 
FSs given in Fig. 1(a). All results presented in Fig. 2 scales 
linearly with the SC gap magnitude as the relationship given 



(al) ^l^QpplApl (a2) 



in Eq. 12 is maintained always. We will not include the scat- 
tering matrix-element C [i.e., Eqs. |6] and |7| in this case to 
theoretically study the behavior of all q vectors as a function 
of energy. The computed results follow the similar behavior 
as sketched in Fig. 1, demonstrating that most of the evolution 
of the QPI maps can be understood from the simple energy 
and momentum conservation rules derived in Eqs. [5]- [7] 

All the intraband scattering vectors lie so close to the strong 
elastic peak at q = that it is often difficult to distinguish. We 
have taken broadening to be 77 = 1 meV which is sufficient for 
the quasiparticle states at = | | to extend up to = | |, 
allowing some elastic scattering at q2 in Fig. 2(al) (although 
it is prohibited in a clean limit). But our small broadening 
does not create visible intensity at all other interband vectors, 
although the experimental data at the corresponding energy 
shows some finite intensity at them, compare Fig. 2(al) with 
corresponding experimental data in Fig. 2(bl). 

At ^^Qpi = |Aq,| ^ |A^|, all interband scattering vec- 
tors between a to 7 FS appear on the QPI map as shown in 
Fig. 2(a2). As at ^Qpi, q2 vector can only show up at ^Qpi 
due to residual broadening of the /3 states upto 7 pocket. The 
separation between qi and q2 can be studied more clearly in 
the overdoped region of hole doped side where the areas of 
the a and (3 bands are distinguishably different. 

The experimental data in Fig. 2(b2) also shows all the cal- 
culated q vectors. Subtle discrepancies in the relative inten- 
sities of each q vector is expected, because we have not in- 
cluded any matrix-element M, and the scattering coherence 
factor C in this calculation. Furthermore, the magnitude of 
each q vectors do not match quantitatively with our calcula- 
tion as the calculation is done at a different doping than the 
experimental data, although in both cases all three FS pockets 
are present. As both QPI and INS data are available as a func- 
tion of momentum and energy for 1 1 1 family LiFeAs, we also 
compare our theory with these data. As shown in Fig. 2(cl) 
and 2(c2), the experimental data"^^ qualitatively agrees with 
the QPI vectors presented in the top panel. In performing 
such comparison, we have to pay attention to the difference 
in FS topology and SC gap symmetries between 122 and 111 
families. 
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FIG. 2. (al) Computed QPI map at Qqpj = |A/3|. As sketched in 
Fig. l(bl), at this energy the q vectors come from intraband scatter- 
ing within /3— FS which concentrate near ^ = 0. The weak inten- 
sity at q2 comes from finite broadening of the quasiparticle states at 
7-pocket which allows some spectral weight of this band to appear 
around ^^qpi = |A^|. (a2) Same as (al) but at Hqpi = |Ac| ~ 
I A^ I . The results agree qualitatively with an earlier one,^^ even after 
neglecting the orbital matrix element term in our case. However, we 
do not see the QPI pattern to be present apart from these two energy 
scales, unless a large lifetime broadening is used in the calculation 
to smear out the quasiparticle states. (bl)-(b2) Corresponding ex- 
perimental QPI data at these two energy scales for electron doped 
pnictide at x = 0.06 taken from Ref. 25 Experimental data show 
all four q vectors depicted in Fig. 1 at two energies, although their 
relative intensities are representative of our calculations in (al) and 
(a2), respectively, due to low experimental resolution. As our phe- 
nomenological calculation do not capture the actual intensity of QPI 
maps, we normalize all QPI maps to their maximum. (cl)-(c2) Ex- 
perimental data^^ for a 111 compound LiFeAs at two bias energies 
^QPi =-1.9 meV in (cl), and Qqpi =-11.7 meV in (c2). Despise 
the differences in the FS topology for this 111 and our 122 system, 
we see that both results are in qualitative agreement. The experimen- 
tal value of A for LiFeAs is 10 meV, which makes the result (c2) 
slightly above the SC region.-^ 



Once magnetic field is applied, the relative intensity of qi 
with respect to that of q3^4 evolves with the strength of the 
field, it has been observed in Fe(Se,Te) compounds, in accord 
with our calculations.^^ 

The proximity of q4 to the reciprocal vector (27r, 0) and 
its equivalent directions, which has been observed both in 
pnictide^^ as well as in iron-chalcogenide^^, has been ar- 
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(a) ^^Ns^lApl +IA,I (b) ^Vs=IAal +IA,I 




FIG. 3. (a) Computed INS map at QJns = | A/3I + | in 1 Fe unit 
cell notation. The intensity peak shifts away from (tt, 0), implying 
that the corresponding resonance peak is incommensurate due to the 
shape of the FSs. (b) Same as (a) but at ^^Jns — |Aa| + |A^|. qi 
peak is closer to the commensurate vector than q2 in (a) as a-FS 
pocket is smaller than the /3-one. The q values are strongly depen- 
dent on the FS areas, hence on doping. The results are consistent with 
the calculated spin-excitation dispersion plot presented in Ref. 12 in 
2 Fe unit cell. Both the incommensurate and the commensurate reso- 
nance peaks are observed in INS experiments.^^ Note that the QPI 
maps shown in Figs. 2(al) and 2(a2) will match exactly with (a) and 
(b), respectively, if a strong magnetic field is applied in the former 
case to eliminate the scattering of same sign of the SC gaps, (c)- 
(d) Experimental data for LiFeAs is presented at Q}jvs=5 meV and 
njj^g=\0 meV from Ref. 44 The data is rotated by 45"" as it is 
avaiable in 2 Fe unit cell notation, while the present calculation is 
performed in IFe unit cell notation. 



gued theoretically to arise from Bragg peak, instead of QPI 
scattering."^^ Comparing the evolution of the intensity of q4 at 
two energies ^^Qpi in Fig. 2(al) and ^gpi in Fig. 2(a2), we 
can deduce that it is not associated with Bragg peak as the 
latter does not have any energy dependence while q4 has. 



B. Computed INS maps of iron-pnictide 

The INS spectra, a direct measure of is calculated using 
Eq.[9]and the results are shown in Figs. 3(a) and 3(b) at two 
energies where INS spectra is finite. The present phenomeno- 
logical approach does not include the overlap matrix-element 
M and the RPA correction. We have shown earlier in Ref. 12] 
that they do not change the essential features of the INS for 
these systems, of course except intensity. 

q2 and qi appear at ^i^g and ^i^g, respectively as ex- 



pected, because they involve sign change of the SC gaps be- 
tween the initial and final states. q2 being smaller than qi, 
will lead to an incommensurate resonance while the latter 
is close to the commensurate one (both are doping depen- 
dence as described above). Both the commensurate"^^ and 
the incommensurat JSI resonances have been detected by INS 
measurements, although the simultaneous presence of the two 
modes is yet to be detected in future measurements with better 
experimental resolution. The experimental data^^ for LiFeAs 
in Figs. 3(c) and 3(d) is consistent with the energy dependence 
of the q2 and qi vectors. 

As mentioned earlier, the INS maps will correspond to the 
QPI maps if the latter is performed at zero magnetic field. For 
the experimental data of LiFeAs, such a compasion between 
Figs. 2(cl) and 2(c2) for QPI patterns and Figs. 3(c)-(d) for 
INS reveals good agreement. 



IV. LAYERED IRON-SELENIDE 

FS properties: -We turn next to the layered AyFe2-x^^2 
based systems. These materials host only electron pockets 
at M point with no hole pocket.^ To model such electronic 
structure, we employ two band tight-binding calculation of 
the t2g orbitals of Fe d bands and the parameters are obtained 
by fitting to the material specific LDA dispersion.^ The result- 
ing FSs are shown by green lines in Fig. 4(a), which match 
well with ARPES FS.^^ 

SC gap properties .--The absence of hole pocket results in 
strong nesting between the electron pockets along qa. We 
have shown earlier that such nesting lead to nodeless and 
isotropic dx2_y2 — pairing symmetry [in 2 Fe unit cell the pair- 
ing symmetry becomes dxy^] and a spin resonance near the 
commensurate vector qa.^ The spin-resonance has recently 
been found experimentally in this class of material by Park 
etaim. 

QPI and INS spectra: -Due to one SC gap and one FS (con- 
centric electron pockets), both the QPI and INS maps appear 
only at one energy scale. The QPI map in this material as 
shown in Fig. 4(bl) resembles the QPI map for pnictide in 
Fig. 2(b2), with the exception that the electron-hole scatter- 
ing q2 is absent here. On the other hand, the INS spectra 
[Fig. 4(c)] is rotated by 45^ in comparison with the same plot 
of pnictide in Fig. 3(b) as the 'hot-spot' for sign-reversal SC 
gap is now aligned along qs. It is interesting to note that al- 
though qs scattering is also present in pnictide as seen in QPI 
maps, but the leading nesting shifts along qi, giving a very 
different pairing symmetry. This observation leads us to con- 
clude that subtle differences between the nesting and the scat- 
tering which can carry rich physical insight can be untangled 
more clearly by comparing the QPI with INS spectra. 



A. Magnetic field dependent QPI maps 

As mentioned before, another way to distinguish the scat- 
tering and nesting vector in QPI map is to study its mag- 
netic field dependence of it at fixed energy. At zero magnetic 
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FIG. 4. (a) Computed FS for layered iron-selenide KccFe2Se2 system 
within two band tight-binding model at a representative in 1 Fe unit 
cell notation. The black to white background depicts the d^2_y2- 
pairing symmetry which takes the form of cos (kxa) — cos (kya) in 
1 Fe unit cell and dxy-wave in 2 Fe unit cell (not shown).^ ^ The ar- 
rows give two interband scattering channels that survive in this class 
of material, compared to iron-pnictide in Fig. 1(a). (bl) Sketch of the 
QPI map according to Eq.[5](the associated scattering matrix-element 
which determines the sign of the SC state of the initial and final states 
is ignored here). (b2) Computed QPI map shows all q- vectors pre- 
dicted in (bl). The small differences in the magnitude of q- vectors 
from the two concentric electron pocket is not distinguishable due to 
finite broadening, (cl) In the INS spectra, q4 does not show up as it 
connects the quasiparticle state of same sign of SC gap. (c2) Corre- 
sponding computational result of INS spectra at QiNS = 2Qqpi. The 
resonance spectra is commensurate at (tt, tt), while the inclusion of 
orbital matrix-element M shifts it to a slightly incommensurate one 
at a critical value of the interaction [/.^ Our prediction agrees with 
recent INS measurements in this class of materials.-^ 



field, the scattering matrix element C only allows scattering 
of quasiparticle states of opposite phase of the SC gap, see 
Eq. |6] Therefore, the resulting QPI map at I^qpi will resem- 
bles the INS spectra at I^ins = 21]qpi, although the q = 
elastic peak will be difficult to remove in the former case. 

At finite magnetic field, the other scattering channel, q4, at 
which the scattering from states of same sign of the SC gap be- 



(a) B-0 (b) B=Finite (c) (a)-(b) 
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FIG. 5. (a) Computed result of QPI pattern in Ka^Fe2Se2 using Eqs.|5] 
and [6] which mimics the zero magnetic field (B) condition at which 
those QPI vectors become illuminated which correspond to the scat- 
tering of states having opposite sign of the SC gaps.^^ (b) Same as 

(a) but using Eqs.|5]and[7] This means, scatterings of same sign of 
the quasiparticle states are only included which mimic a finite ap- 
plied magnetic field condition. In principle, a magnetic impurity has 
a scalar component which will allow some scattering for qs vector as 
well but relatively weak in intensity, (c) Differences between (a) and 

(b) . The positions of maxima and minima corresponds to qs and q2, 
respectively, an unambiguous method to find out which scattering 
channel involve sign change of the SC gap, a test which will confirm 
the presence of d-wave gap in these materials. 



comes turned on as shown in Fig. 5(b). In practice, a magnetic 
impurity is always associated with a scalar component (not 
considered in the present calculation for simplicity), therefore, 
qs should also be present. The relative intensity can be moni- 
tored by tuning the strength of the magnetic field. On the 
otherhand, the difference between QPI maps in Fig. 5(a) and 
5(b), as shown in Fig. 5(c), can also be used to identify the 
sign-changing 'hot-spot' vector. The positions of maxima and 
minima corresponds to qs and q4, respectively. The experi- 
mental detection of the maxima and minima can be tested to 
confirm the presence of d-wave gap in these materials. 



V. CONCLUSIONS 

We have developed a phenomenological approach that al- 
lows one to connect the INS and QPI maps. We point out that 
both the INS and QPI maps arise due to the inelastic and elas- 
tic scattering of the Bogoliubov quasiparticles, respectively, 
and bear a one-to-one correspondence as a function of scatter- 
ing vector and energy. In fact, inelastic and elastic scattering 
complements these two spectroscopies to quantitatively and 
unambiguously identify the nature of pairing symmetry in un- 
conventional superconductors. By applying our model in iron- 
pnictide and layered iron selenide compounds, we show that 
the QPI maps at zero magnetic field corresponds exactly to 
the INS spectra at their representative energy (the q = elas- 
tic scattering and other spurious effects in QPI maps can be 
ignored because of "contamination from Bragg peaks"). We 
also show that upon applying magnetic field , the QPI scat- 
tering at same sign of the SC gaps can be illuminated in lay- 
ered iron-selenides and produce detectable changes in QPI. 
We point that evolution of the QPI maps can be implemented 
experimentally to test the possible nodeless d-wave pairing in 
this class of materials. 
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